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Abstract 

An analysis of the real and imaginary part of the polarization po- 
tential in terms of the relative contributions of the single collective 
states for the 20S Pb -)- 208 Pb system has been done. The polarization 
potential has been calculated within the Feshbach formalism taking 
into account the collective states calculated with the Wigner function 
moments method. The contribution of the isoscalar giant dipole reso- 
nance states has been estimated being of the order of 10 - 20 % of the 
total at relatively low incident energy. 



PACS : 25.70. G, 25.70. C 
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INTRODUCTION 



In the collisions between complex nuclei, the excitation of the internal degrees 
of freedom gives rise to a modification of the bare nucleus-nucleus potential, 
calculated for instance by the double folding procedure This modification 
is usually referred to as polarization potential. Many authors ]2j, ||, £|, §], 
using different approaches and various prescription to calculate the form 
factors, have studied the relevance of the excitation of vibrational collective 
states (both low lying and giant resonances) for the determination of the 
polarization potential. The role played by the nucleon transfer process has 
also been investigated |J. 

Very recently, by using (a, a'7) reactions, some experimental evidence 
has been found || on the existence of low lying isoscalar dipole modes in 
208 Pb, 90 Zr, 58 Ni and A0 Ca. These modes are of compressional nature and 
are quite well described by the method of the Wigner Function Moments 
(WFM) J7| . Due to their isoscalar nature, it is interesting to investigate how 
much they contribute to the polarization potential. In order to calculate the 
latter within the same microscopic approach used in previous studies || |J 
one needs the transition densities of the considered excited states. As shown 
in the next section, these can be extracted from the WFM results by using the 
linear response theory. Calculations performed for the 208 Pb + 208 Pb system 
at various incident energies show that the contribution of these isoscalar 
compressional modes to the polarization potential amounts to 10-20 %. 



2 DESCRIPTION OF WFM METHOD 



The starting point for this method is the time-dependent Hartree-Fock equa- 
tion for the density matrix p g (ri,r 2 ,t): 



index for protons (q = p) or neutrons (q = n). A Skyrme-type effective force 
(SGII) is used as nucleon-nucleon interaction. 

Equation (1) is transformed to an equation for the Wigner function 
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^ = fsm[|(Vf.V;-Vf.V;)]^ ) / 9 , (2) 

where H$(r,p) = J e - ip s / n (r + s/2 | | r - s/2) ds and r = (n + 

r 2 )/2, s = ri — r 2 . The moments of f q in momentum space define the 
nucleon densities n q (r, t), the mean velocities u q (r, t) and tensors of any rank 

Pqij ...fc(r 'j t) ■ 

n q (r,t) = J /,(r,p,t) dp, 

m-n g -Ug(r,t) = J p /,(r, p, t) dp, 

Pqh...in = ml ~ n J (Ph - mu qil ) . . . (p in - mu qin )f q (r, p, t) dp 

where m is the nucleon mass. The method consists in taking moments of 
equation (2) in phase space, i.e. in integrating Eq.(2) in p and r spaces 
with different weights, in order to obtain a set of dynamic equations (virial 
theorems) for some chosen moments of the Wigner function. For example, 
to describe the negative parity states 1~, 2 _ , 3~ one must integrate Eq.(2) in 
momentum space with the weights l,Pi,PiPj,PiPjPk (Pi are components of the 
momentum p) as a first step. One arrives at a set of coupled time evolution 
equations for the density (continuity equation), for the mean velocity u q i(r, t) 
(Euler equation), for the pressure tensor P qi j{v,t) and for the third rank 
tensor P qi jk(r,t). Then the obtained equations, weighted with X{ and XiXjXk 
for n q , with 1 and XjXk for u q i, with x q for P qi j and lastly with 1 for P q ijk, 
are integrated in r-space. Thus one arrives at a closed system of coupled 
dynamical equations for different integral characteristics of the nucleus - 
Cartesian tensors of the first and third ranks: / n q (r, t) Xi dr, J u qi (r, t) dr, 
J n q (r,t) XiXjXk dr, J u qi (r,t) XjX k dr, J P qij (r,t) x k dr, J P qijk (r,t) dr. 

In the case of simple interaction (a harmonic oscillator with a separa- 
ble multipole-multipole residual interaction, for example) the integrals con- 
taining the interaction can be expressed in terms of the Cartesian tensors 
mentioned above. As a result one obtains a system of nonlinear equations 
for the integral characteristics of different multipolarities. In the general 
case of a realistic interaction the integrals containing the interaction can 
not be expressed in terms of integral characteristics without any approx- 
imations. This problem is solved rather easily in the case of small am- 
plitude motion. To obtain the corresponding equations we vary the virial 



3 



equations taking n q (r,t) = n q °'(r) + dn q (r,t), u q i(r,t) = u qi (r) + 5u q i(r,t), 

P qij (r,t) = J$](r) + SP qij (r,t), P qijh {r,t) = pj® k (r) + SP qijk (r,t) and ig- 
noring terms quadratic in the variations 5. So we arrive at a set of linearised 
equations, which are a very convenient means to study the collective small 
amplitude motion. These equations are simple differential equations. In the 
case of harmonic oscillations they become algebraic equations. The coeffi- 
cients in these equations depend only on the ground state properties. They 
are linear combinations of integrals over the nuclear volume of different pow- 
ers of ground state particle densities n q °\ kinetic energy densities P^l and 
of their space derivatives. 



3 TRANSITION DENSITY 

The transition density is defined as the matrix element (0|p(r)|a) of the 
density operator 

N 

p(T) = J2H*-h) (3) 

i=l 

where is the position operator of particle i and |0) and \a) are the ground 
and excited state, respectively ||. The problem of WFM method consists in 
finding this matrix element without calculating the wavefunctions of |0) and 
| a). This can be solved with the help of linear response theory. The linear 
response to the perturbation operator 0{t) = Fexp (—iujt) + i^exp (+iut) 
can be written as 



Pki 



(i) _ I x ^ (- (0\aja k \u)(u\ala q \0) (0\ala q \i/)(i/\a}a k \0) ^ 
(0\aja k \u)(u\F\0) (0\F\v) (u\a\a k \0) x 



UJ 



Vty UJ + Fly ' 



where a\,a k are creation and annihilation operators, Vt u the eigenfrequen- 
cies of the system, and F = Y^, pq fpq^ p a q . Using the second quantization 
representation for the density operator 

P( r ) = Y. d pi( v ) a l a n ( 5 ) 

qp 
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where d pq (r) = (p\5(r — r)\q) = 0*(r)0 g (r) and (f) q (r) being a basis of single 
particle wavefunctions we can write the following equation for the change of 
density 



Sn(r) = J2pk } di k (r) 
ki 

1 ( (0\p(r)\u)(u\F\0) (0\F\u)(u\p(r)\0) } 

which we calculate with the help of the WFM. By varying the continuity 
equation 

— = -div {n q u q + r]n q n q/ (u q - u' q )} (7) 
we express the change of density in terms of small displacements £, 



qi 



3 Q 

6n g = ~ X, { n g£qs + Vn q n q ,(£ qs - ^ s )} (8) 

where ^ q i(r,t) = 5u q i(r,t). These displacements are parametrized in the 
following way 

3 3 
£ qi (T, t) = L qi (t) + L qi,A t ) X i + L Qi,jk(t) XjX k (9) 

j=l j,k=l 

The tensors L qi (t), L qi j{t), L qi jk(t) are found by solving the system of cou- 
pled dynamical equations involving moments of the Wigner function. Their 
derivation was described in ref. |7], M. It is evident that one can extract the 
information about the transition density from eq. @ taking the limit 

lim h(u - tt a )5n(v) = (0|p(r)|a)(a|F|0) . (10) 

U) ►Ha 

Thus, to determine the transition density (0|p(r)|a) from this expression it 
is necessary to know the matrix element (a|F|0). By using eqs. (P) and (|10|) 
one can get the square modulus of the required matrix element : 

lim ft(w-fi a ) [ dr/(r)£n(r) = (0\F\a}(a\F\0) = |(a|F|0)| 2 . (11) 



5 



which is true for any hermitian single particle operator F. In the case of a 
real operator one can write 



(a|F|0)=±[ lim h(u - Q a ) [ dr f(r)5n(r)] 1 . (12) 
Then the final expression for the transition density is 

(0|p(r)|a> = . (13) 

[ hm w ^n a ft(a; - ft a ) / dr f(r)Sn(r)\ 2 

Combining eqs. C |X3|) and (||) we get the explicit expressions for the transition 
densities of the various multipolarities 

(0|p(r)|2+) = r ^-(A 1+ A 2 n(r))Y 2 ,, (14) 

(0|p(p)|3"> = r 2 ^-(B 1 + B 2 n(r))Y 3 ,, (15) 
(0|p(r)|l~) = {rn(r)(C 1 + C 4 n(r)) + 

^[C 3 + C 2 r 2 + n(r)(C 5 + C q t 2 )\)Y^ , (16) 
<0|p(r)|0+> = {A(n(r) + ^) + 

D 2 (n 2 (r) + ^rn(r)^)}Foo (17) 

where the constants A and Z) are combinations of the tensors L qi j, B - 
of tensors L qi j k , C - of tensors L gi and L qi j k . These expressions are too 
complicate, to be written here. So we describe the derivation of formulae for 
the most simple constants (Ai, A 2 ) in the Appendix, just to give a general 
idea of the procedure. 

The transition density of Giant Octupole Resonance (GOR) differs from 
that of Low Energy Octupole Resonance (LEOR) by the values of the con- 
stants B. It is necessary to note that in our calculations LEOR includes also 
the contribution from the lowest 3~ state, because the WFM method gives 
the centroid of all 3~ states lying below GOR. The expressions ( fl4T) , fllBl) for 
the transition densities of 2 + and 3~ excitations differ from that of Tassie 
model by terms quadratic in the density (terms in A 2 and B 2 ). Numerically 
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these terms turn out to be small. The formula for the transition density of + 
practically coincides with the one derived in ||T^]. The only difference is the 
term proportional to D 2 , which however is numerically very small. The ex- 
pression (16) can be compared with formula of Harakeh and Dieperink [|TT|1 , 
which was derived supposing that all energy weighted sum rule (EWSR) is 
exhausted by one state. All the three terms linear in the density have the 
same radial dependence, but the corresponding coefficients Ci, C 2 , C 3 have 
nothing in common with that of ref. |TTJ. Moreover, the terms proportional 



to n(r) 2 give a noticeable contribution in this case. So the total radial depen- 



dence of 1~ transition density may be different from that of |TT|. We have 
three low- lying 1~ states contributing to the isoscalar EWSR. Their proton 
and neutron transition densities are shown in Fig.l, left and right side, re- 
spectively. The dashed lines give the transition densities calculated by using 
only the linear terms in the density, while the solid lines include also the 
quadratic ones. It is seen that the transition densities of all three excitations 
are rather large inside the nucleus demonstrating their compressional nature. 
The highest excitation has the largest contribution quadratic in the density 
and has a rather big admixture of isovector mode. 



4 POLARIZATION POTENTIAL 

The most elegant way to calculate the polarization potential is to follow the 
Feshbach formalism |2|, In this approach is very easy to take explicitly 
into account the contribution of some definite collective states. Here we will 
briefly describe the method which has been extensively described in ref. [Q, [|]. 
In the Feshbach approach the effective heavy ion interaction is written as 

V(R, R') = (00KR)|00)<5(R-R) 

+ J2 (00\v(R)\K 1 K 2 ) x GWRR'X^i^KROlOO) 

K!K 2 

= V F (R) + AV(R,R') . (18) 

In the first term the nucleon-nucleon interaction v(R) is double folded with 
the ground state densities of the two nuclei. In the second term the sum is 
over all the non-elastic channels, and v (R) is double folded with the transition 
densities of the two nuclei: it describes the coupling of the elastic channel 
to the non-elastic ones. This term is the so called dynamical polarization 
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potential. Its physical meaning is transparent from eq. flIB|) : The interaction 
acting at the distance R' takes the system in one of the eliminated channels, 
then it is propagated at another distance R where the interaction, acting 
again, brings back the system into the elastic channel. Then the polarization 
potential is non-local, and if one or more of the eliminated channels is open 
it is also complex and its absorptive part describes the loss of flux from the 
elastic channel. The couplings make AV also energy dependent because of 
the appearance of the energy in the propagator Gk 1 k 2 (^ R')- 

The calculation of eq. (|Tj|) is a very difficult task because of the presence 
of the propagator, then we calculate it in the WKB approximation: 

n ( nB \~ E exp(iM KlK2 (p)s) 

G Ki k 2 {P, s) ~ -7TT2 ( 19 ) 

2irn s 



with 

M 2 KlK2 (p) = p [E cm - E Kl - E K2 - V L (p) - V c {p)\ (20) 



where 

p=^(R + R'); s = |R-R'| (21) 

The local optical potential V^(p) is given by 

V L (p) = V F (p) + AV L (p) (22) 

while Vc is the Coulomb potential between the two nuclei. A procedure 
to define the local polarization potential AV^ will be given later on in this 
paper. 

Within the double folding approach the form factors can be written as 

F^o(R) = (00KR)|ATi0> = J dndr^^r^v^ - r 2 + R|)p 0Q (r 2 ) (23) 

where p Ki0 and p 00 are the transition density and the ground state density, 
respectively. These are calculated within the WFM method, as described in 
section |^. The details of the calculation of AV(p, s) is reported in ref. 0. 

The potential so obtained is non-local, but since the range of non-locality 
of AV is small with respect to its radius [|J we can use a standard proce- 
dure |H| to obtain a local potential from a non-local one. Then 

AV(p) = 4tt J Jo (ks)ReAV(p,s)s 2 ds (24) 

W(p) = 4vr J Jo (ks)ImAV(p,s)s 2 ds (25) 
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where j is the Bessel function and 



2/i 

h 



2 [E cm - V F (p) - V c (p)} . (26) 



5 RESULTS AND DISCUSSION 

We have done calculations for the system 208 Pb + 208 Pb at several incident 
energies. The levels used in the calculations are reported in Table I : they 
were obtained with the WFM method. We have parametrized the ground 
state density with a Fermi distribution n(r) = no/(l + exp [(r — R)/a\) whose 
parameters are the following [13]: R = 1.115 A 1 / 3 — 0.53 A^ 1 ^ 3 fm, a = 



0.568 fm, while no is fixed by the condition 4-7T / °° n(r)r 2 dr = A. As already 
stated in section |3], the transition densities have the advantage to have an 
analytical form from which one can see that their first term corresponds to 
the Tassie model in the case of 2 + , 3~ and + multipolarities. The corrections 
with respect to the Tassie ones are not important. Formula for 1~ is new and 
has nothing common with Tassie model, because it takes into account the 
center of mass movement. The bare potential has been obtained by double 
folding the nucleon - nucleon effective interaction M3Y with the ground 
state density of the two nuclei. In the same way we have obtained the form 
factors eq. ( [23] ) by double folding the M3Y with the ground state density of 
one nucleus and the transition density of the other one. Examples of form 
factors are shown in Fig. 2: on the right side we show the one related to 
the LEOR, while on the left side are reported the ones concerning the three 
isoscalar GDR states. In the latter case the oscillations are due to the fact 
that the form factors change sign around 10/m. 

In order to show the relative contribution of the states included in the 
calculations we have computed the polarization potential for various energies 
as function of the relative distance R. In fig. 3 we show the contribution due 
to different multipolarities( as indicated in the figure) to the real (left) and 
imaginary (right) part of the polarization potential plotted in percentage of 
the total potential as function of the relative distance R. We notice that the 
most important contribution to both real and imaginary part is due to the low 
lying 3~ state. As the incident energy is increasing the relative contribution 
of the GQR becomes more important. This different behaviour has already 
been found in previous studies || ffl. The novel result is represented by the 
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contribution of the isoscalar GDR states which lies between 10-20 %, that is 
to say comparable to the contribution of the GQR at least for the imaginary. 

We have also calculated the polarization potential at a fixed relative dis- 
tance (R = 16 fm) as a function of the incident energy taking into account 
the explicit contribution of the different multipolarities. The results are 
shown in fig. 4 where the different curves correspond to the contribution of 
the different states as indicated in the figure. Again we notice the importance 
of the LEOR, which is overwhelming with respect to the others states in a 
large range of incident energy for the case of the imaginary part, while for 
the real part the contribution of the GQR, at higher energies, is essentially 
the same than the one of the low lying 3~ state. The contribution of the 
isoscalar GDR to both real and imaginary part of the polarization potential 
is more important (at this distance) at low energy, and it decreases with the 
energy slower for the absorbitive part than for the real part. 

The interplay between the LEOR and the GQR states seems to contradict 
the findings of ref . ||, . In particular the contribution of the LEOR to the 
absorbitive part is very big with respect to the one of the GQR. As explicitly 
stated in the second paper of ref. [Q, this is due to the fact that we are 
using different transition densities (hence different form factor) from the one 
used in ref. Q which were calculated in the framework of the Random Phase 
Appr oximat ion . 

6 CONCLUSION 

We have done an analysis of the real and imaginary part of the polarization 
potential in terms of the relative contributions of the single collective states 
for the system 208 Pb + 208 Pb. The polarization potential has been calculated 
within the Feshbach formalism taking into account the collective states cal- 
culated with the WFM method. Within this method it is possible to obtain 
also isoscalar GDR states whose contribution to the polarization potential 
has been estimated being of the order of 10-20 % of the total at relatively 
low incident energy. Then this contribution should in principle show up in 
some experimental observable, like for instance the elastic cross section. The 
other multipolarities give a contribution which has been analyzed in great 
details in previous work 0, || where slightly different results were obtained. 
This can be ascribed to the different form factor we have used in the two 
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cases. 
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APPENDIX 
Transition density for 2 + excitation. 

Multiplying (8) by r 2 y 2/J (6 l , ip) and integrating over angles we find the 
radial dependence for the change of density 

Ati Oti 3 
Sn q (r) 2fl = -J^{ r ~Q^ ■ £q2i* + W-Q^{n q ■ ry) • (£ ?2jU - , {Al) 

where the components of the irreducible tensors C q2ll are linear combinations 
of L q ij. In the case of a spherical nucleus we can study any component of 
£> q 2n, for example C q2 2- 

£q22 = L q ii — L q 22 + 2iL g i2, L^j = Lij + Lj t i. (A2) 

For a spherical nucleus L q n = L q2 2 so we will concentrate our attention 
on L q i2- Following the prescriptions of chapter 2 and using formulae (8,9) 
one can derive the system of coupled dynamical equations for L ql2 and 
JSP ql2 dr = n ql2 (see ref. @): 

L q i2 — 2b q \L qV2 + 2b q2 L q 'i 2 — 2b q3 (37r q i 2 = 

-ifi J r ^Y 2 ,-W(t)dr, (A3) 

(37T q 12 + d q i • L q i 2 — d q 2 • L q i 12 = 0, 

where (3 = A ^ q4 and W(t) = e p r 2 (Y 2fl e~ iuJt + Y 2 \t iuji ) is the external 
field. The corresponding system for the other kind of nucleons is obtained by 
changing q <-> q'. The following notations have been introduced here: 



ih A + ->4ki £ * + 2T: " + : ' {1Tl + 2T2 ' ] + 5 " ' 



z t 8 

bq2 = jA + 1-T^Sq + z(7T X + 2T 2q ) + 5q )P ■ —(pTJUoZ 

ma\ An a i 15 
2ixm , 15 ocl 



ft A 8 maf 
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1(1 + 1)^ + 2^) + ^ 
+2a(a -!)(! + 2x 3 )zz'} ■ (111 + 2^ +1 ) , 



A = j(l + ^)(St + 2 V nl) + - [ ^{(a + l)Ml-x 3 ) + 2(2 + x 3 )] + 



£ q = z(KT + 2 V Kl /3 ), v = 7 rn i^, 

m A 

,dn 



a»= n[rYr v dr, Vf v = / n(r)V(— f dr 

•J "0 C/f 

Jo or Jo o z r 



n(r) =n p (r) +n n (r), z =j£, Z ' = ~A' V= ~W' f+ = *i + *2, 

i_ = ti — t 2 . The quantities t 0) ^2, £3? cr, ^0; ^3 are the parameters 
of Skyrme forces, e p the proton charge, m the nucleon mass. Deriving these 
formulae for b q d q iwe used the Thomas- Fermi approximation for the ground 
state pressure tensor P qi j(r) = 5jj7?i ? (r) 5 / 3 with 7 = J^(37r 2 ) 2//3 . 

The system of differential equations (A3) become algebraic when one 
takes into account the evident time dependence of variables L ql2 (t) = L q \ 2 e iwt 
and n ql2 (t) = ii q i 2 e llAlt . the solution of this system is : 

L ql2 = A q /A, L q , 12 = A g ,/A, (AA) 

where A is the determinant of system (A3). The determinant A q is obtained 
by substituting the corresponding column of A by the right hand side column 
of (A3). A is a biquadratic polynomial and has two roots Q 2 and f2 2 which are 
interpreted as isoscalar and isovector giant quadrupole resonances. Writing 
the determinant in the form A — c - {oj 2 — VL 2 ){uj 2 — Vt 2 ), where c is a constant 
and putting (A4) into (Al) one can calculate easily the limits in (13) and 
find the constants Ai, A 2 . 
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FIGURE CAPTION 



Fig.l 

Proton (left) and neutron (right) transition densities for low-lying dipole ex- 
citations in 208 Pb with the energies E x = 7.65MeV, E 2 = 8.99MeV and 
£3 = lO.OlMeV , from top to bottom. The dashed curves are the results of 
calculations with only the linear density terms (formulae (14)- (17)). 

Fig.2 

Form factors for 208 Pb for the reaction 208 Pb + 2 08p^ f or ^ wo different mul- 
tipolarities. The three curves in the left part correspond to the isoscalar 
GDR states with excitations energies E x = 7.65 MeV (solid curve), E 2 = 
8.99 MeV (dashed curve) and E 3 = 10.01 MeV (long - short dashed curve). 

Fig.3 

Contribution of the different multipolarities ( as indicated in the figure) to 
the real (left) and imaginary (right) part of the local polarization potential 
plotted in percentage of the total potential as function of the relative distance 
R for three different values of incident energies for the reaction 208 Pb+ 208 Pb. 

Fig.4 

Contribution of the different multipolarities ( as indicated in the figure) to 
the real (upper part) and imaginary (lower part) part of the local polarization 
potential plotted as function of incident energy for a fixed relative distance 
R = 16/m for the reaction 208 Pb + 208 Pb. 
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Table 1. 

Isoscalar collective states of 208 Pb used in the calculations. 



J" 


E, MeV 


EWSR, % 


\- 


7.65 


5.79 




8.99 


2.94 




10.01 


1.13 


2+ 


11.82 


33.35 


3" 


4.65 


15.14 




22.21 


16.20 


0+ 


13.63 


34.28 
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